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Abstract A technique to find Pareto-optimal solutions to multiple objective linear programming problems under 
intuitionistic fuzzy environment is presented in this paper. In 1997, Angelov proposed optimization technique 
under intuitionistic fuzzy environment. In 2009, Jimenez and Bilbao showed that a fuzzy efficient solution may 
not be Pareto-optimal solution in case that one of the fuzzy goals is fully achieved. In 2015, Wu, Liu and Lur 
redefined membership function of fuzzy set theory and proposed another two phase technique to find Pareto- 
optimal solution. In this paper it is observed that strictly monotonic part of both membership and non- 
membership functions of intuitionistic fuzzy goals are useful in optimization technique. Further, better optimal 
solutions may be found if constraints that sum of minimal level of acceptance and maximal level of rejection not 
exceeding unity as well as minimal level of acceptance not exceeding unity are removed. Moreover few such 
constraints, used in existing techniques, may make a problem infeasible. Consequently, new functions: T (+) - 
characteristic functions and ^'-characteristic functions are defined in place of membership function and non- 
membership function respectively in intuitionistic fuzzy decision making and used in proposed algorithm. The 
concept of Pareto-optimality in intuitionistic fuzzy environment and new method to test Pareto-optimality of 
solution is proposed in this paper. Necessary counter examples are given and one application on production 
industry further illustrates proposed algorithm. Conclusions are drawn at last. 

Keywords Intuitionistic fuzzy optimization • I-Pareto-optimal solution • Multi-objective optimization • Pareto- 
optimal solution • Production planning • T (+) -characteristic function • T ( ' -characteristic function. 

1 Introduction 

In literature, the framework of multiple objective linear programming problems (MOLPP) with k objective 
functions Z ■ (x) = c-x, i = l ... k is usually considered as the following model 

T T 

Min z(x) = (z^(x), Z2 (t),...,z^.(x)) -(c^x, c^x, (l) 
subject to x e X 
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where X ={ig R" : Ax <b,x>0},b = (b l ,b 2 ,...,b m )e R m and A is a mxn matrix and R denotes the set of 
real numbers. By assuming that the decision maker (DM) has imprecise aspiration levels for each of the 
objective functions Z ; (x) = C t X , i = 1. . .k of model (1), several methods have been proposed in literature to 
characterize Pareto-optimal solutions to MOLPP (1). 

In one such approach, fuzzy set theory is used [see. Bellman and Zadeh 1970]. Fuzzy optimization technique 
is more flexible than deterministic optimization technique and allows us to find solutions that are more adequate 
to real life. And one popular method for finding Pareto-optimal solution under fuzzy environment is two phase 
method. Here max-min operator is usually applied due to its easy computation but it is not guaranteed to yield 
Pareto-optimal solution if more than one optimal solution exists. Using the linear membership functions of 
fuzzy objective functions (and constraints till they are symmetric), max-min operator proposed by Zimmermann 
(1978) can be applied to solve model (1) and find the maximum membership grade that each of the membership 
functions can simultaneously attain. Thus model (1) can be converted into crisp linear programming problem as 
follows [see, e.g. Tanaka et al 1974, H. Zimmermann 1985] 

max A (2) 

subject to 

Mi(Zi(x))>Z, i = 1 . . . k, 

0<A<1, xe X. 

There exist several methods in literature to obtain optimal solution to this problem (2) [see, Guu and Wu 
1999]. Sakawa (1987) contained a rather comprehensive survey regarding to interactive methods for MOLPP 
and fuzzy MOLPP [see Sakawa 2013], Dubois and Fortemps (1999) proposed a multi-step procedure instead of 
two phase method to solve model (2). Recently Jimenez and Bilbao (2009) suggested that fuzzy-efficient 
solution may not be Pareto-optimal solution in case that one of fuzzy goals is fully achieved. And their proposed 
procedure extended two-phase approach of Guu, Wu (1997, 1999) and approach of Dubois, Fortemps (1999) to 
attain Pareto-optimal solution to MOLPP under fuzzy environment. But according to Yan-Kuen Wu et al 
(2015), proposed approach by Jimenez and Bilbao (2009) cannot guarantee to be one general procedure to attain 
Pareto-optimal solution to model (1). 

Yan-Kuen Wu et al (2015) redefined membership function of fuzzy set theory. But that definition of redefined 
membership function jU- (c-x)\fi = \...k and its usage in mathematical model as well as given numerical 
examples in their paper are not analogous. It is clear that only part of redefined membership function 
jU- (c-xj\/i = \ ...k where value of objective function does not exceed sum of goal and tolerance (for 

minimization type of objective functions) is used in mathematical model and in given numerical examples. But 
these were not added to set of existing constraints. Thus in problem formulation, redefined membership function 

ft- (c-x)\/i = l...kis strictly monotone where as it is not so in definition. Only upper bound of membership 

function was removed but there still exists lower bound of membership function at zero [see, Wu et al 2015], As 
already observed by us, in many cases, if the lower bound of membership function remains intact, the problem 
may become infeasible. 

Deep analysis by us further highlights the cause to lie in definition of membership function of fuzzy set itself. 
Membership function of fuzzy set provides satisficing result when extreme ends of imprecise information can be 
quantified within boundary of zero and one. But it may not be always logically correct to quantify or measure 
each kind of imprecise information within one bounded subset of real line. 

Consider fuzzy set A of all countries that may win next Soccer World Cup in 2018. Most of us agree with the 
prediction that Japan will not be the champion. So it can be safely stated that Japan is in A with membership 

value 0. Again chance of India’s winning is also nil. So India is also in A with membership value 0. It means 
that Japan and India are at same level! But it is certain that even if Japan wins 2018 Soccer World Cup, 
unfortunately India will not be able to participate and hence she will certainly not be the winner. So membership 
function of fuzzy set theory fails to explain this case! 

An element in fuzzy set may lie partly or never lie or must lie in that set. The idea of lying, lying partly or not 
lying is well quantified / measured by membership function of fuzzy set theory. But it does not suit in case of 
must lying (membership value in fuzzy set is always one), and never lying (membership value in fuzzy set is 
always zero). 
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On the other hand, fuzzy set theory has been widely developed and several new modifications and 
generalizations have appeared. One of them is the concept of intuitionistic fuzzy (IF) sets, introduced by K. T. 
Atanassov (1986). It considers not only degree of membership (z ; (x)), i = l...kbut also degree of non- 
membership V^Zj (x)), i = 1. . . k, such that sum of these values does not exceed unity to each IF objective. So 

MOLPP (1) may be treated as MOLPP under IF environment as well. And analogous to above discussions, in IF 
sets, introduction of two new functions in place of both membership functions and non-membership functions 
are necessary. 

The rest of the paper is organized as follows: Sect. 2 introduces definitions. Sect. 3 highlights loopholes of 
existing technique and suggests remedy to find Pareto optimal solutions to MOLPP under IF environment. Sect. 
4 proposes related algorithm. Sect. 5 presents numerical examples where existing technique fails and proposed 
methods gives acceptable solution. Sect. 5 also presents an application on purchasing bulk raw materials of a 
large-scale integrated steel plant. Sect. 6 gives concluding remarks. 

2 Definitions 

It is observed that an element of an IF set can lie partly or never lie or must lie in that set. Moreover in case of 
decision making under IF environment, one interesting and useful property of membership function is that 
higher value of membership function always gives better result of objective function; and analogous property of 
non-membership function is that lower value of non-membership function gives better result of objective 
function. And usual membership functions, non-membership functions are not the only candidates having those 
amazing characteristics. Hence in IF environment, new functions viz. T (+) -characteristic function in place of 
membership function and ^'’-characteristic function in place of non-membership function are defined. 



Definition 1 

Let S denotes the universal set and a be an IF subset of S. Then T ,+) -characteristic function is denoted by 
7L <+ * (x) and is defined as T~/ +l S — > R that assigns a real number T~/ + ^ (x) to each element xe S , here 

T-} + \x) represents degree of membership or acceptance level of xe S in a . And ^’-characteristic function 



is denoted by T-j" \x) and is defined as T-J' ^ : S — > R that assigns a real number T~j~ 
A 1 A A 1 



(x) to each element 



X e S , here T~ ■ (x) represents degree of non-membership or rejection level of x e S in a . 

A 1 



In model (1) it is usually not possible that all objective functions will simultaneously attain optimal values, 
due to their conflicting nature under IF environment. Concept of optimal solution is thus replaced by I-Pareto- 
optimal solution to MOLPP (1) under IF environment and consequently Pareto-optimal solution. 

Definition 2 

A decision plan X 0 e X is said to be an I-Pareto-optimal solution to the MOLPP (1) under IF environment if 
there does not exist another y e X such that 7j + ( Zj (x 0 )) < 7j + (z ; - (y)), (z ; - (y)) < 7j (z ; - (Xq)), Vi, i ^ j ; and 

Tj + ( Zj U 0 )) < Tj + ( Zj(y )) , T~ (Zj(y)) < T ~ (z . (x 0 )) for at least one j. 



Definition 3 (Interactive fuzzy satisfying method, Sakawa 1987). 

A decision plan X 0 e X is said to be a Pareto-optimal solution to the MOLPP (1) if there does not exist another 
ye Xsuchthat z ; (_y)) < z ; (x 0 )for all i, i^j and z.(y) < Z (x 0 ) for at least one j. 



3 Obtaining I-Pareto-optimal solution and Pareto-optimal solution 

Within the scope of multi-objective decision making theory, Pareto-optimality of solution is a necessary 
condition in order to guarantee the rationality of a decision. Now applying the concept of IF sets; Angelov 
(1997) formulated IF optimization (IFO) problem as follows 
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Fig. 1 Strictly monotonically extended membership function and non-membership function (not in scale) 



max jUj ( Zj (x)), i = 1 . . . k, (3) 

min v,.(z ; (x)), i = l...k, 
subject to 

Mi ( z t (xj) + Vj ( Zt (x)) < 1, i = 1 . . . k, 

(jc)) > v. (z. (x)), i = l...k, 

F,(z,(x))>0,i = 1 . . . k, 
xeX. 



Which was further simplified by Angelov (1997) as follows 

max a- p (4) 

subject to 

Miizpx)) > a, i= 1 . . . k, 
v i {z i (x)) <p,i = 1 • • • k, 
a > p, 
a + p< 1, 

P> 0,ie X. 

where OC denotes the minimal degree of acceptance and f3 denotes the maximal degree of rejection of IF 
objective(s) and constraint! s). Without loss of generality, it may be assumed that all constraints of model (4) are 
active. It is observed that in problem formulations as well as during computations, strictly monotonic parts of 
definitions of membership function and non-membership function are used. And constraints 

M i (z i (x))>aP = l...k,andiA(z,(x))</?,i= l..k, are suitably defined by relations a + P <l,Ct> P , 
P > 0 and hence a > 0 in model (4). It is well known that a triangular IF goal to minimization type of objective 

function is usually denoted by (cv,b,b') where //(zia)) = 1 , v(z(a)) = 0 , jU(z(b)) = 0 and v(z(b )) = 1 with 
b < b as shown in Fig. 1. 

Here the constraint a + P < 1 represents that sum of minimal degree of acceptance and maximal degree of 
rejection cannot exceed unity. But lemma 1 shows that on left side of , strictly monotonic part of definitions 
results in sum of membership and non-membership values to always exceed unity. Hence in presence of 
constraint a + P < 1 , objective function z- (x), i = 1 . . . k, may not attain value less than (l t . But since objective 
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function z ; . (x),i = 1 . . . k, is of minimization type and IF goals and/or tolerances are imprecise in nature and 
any value less than £7 ; to objective function Z- (x),i = l...k, is more acceptable to DM, so presence of 
constraint a + /? < lmay result in less acceptable optimal solutions, as shown by numerical example 5 . 1 . 1 . in 
this paper. Hence it is suggested to remove the constraint a + [5 < 1 so that objective functions Z, (x), i = 1 . . . k 
can attain more preferable values. 



Lemma 1 Let{a i ',b i ,b j ) be triangular IF goal of objective function Z t (jc), i = l...k of minimization type to 

// 

MOLPP (1). Then for any CL < a i , usage of strictly monotonic part of definition of membership and non- 
membership function during computation yields that sum of membership degree and non-membership degree 

exceeds unity i.e. jU i {z j (.Cl i )) + V i (z i (a ! ))>l,Vflt < a r 

Proof Suppose (a i ',b i ,b i )be triangular IF goal to objective function z ; (x),i = l...k, of minimization 
type to MOLPP ( 1) under IF environment. Then V? = 1 . . . k, from definition, it is known that jU i (z ; (flj )) = 1 

and P I (z ! (a ; )) = Oand /l l (z i (b j )) = 0 and V j (z i (b l )) = lwith£> < b. . For all / = l...k , let us choose any 

// rr 

a ; < a ; such that V; = I . . . k,a- = zAx) for some x e X . Then strictly monotonic property of membership 

. * b.-a. * ct;- a" 

function gives fU \Z \Cl )) = and non-membership function gives V {z (Cl- )) = . 

b t ~a t b'-a . 

Now V/ = 1. . . k, we have 



M i U i (a")) + V i (z i (a'')) 
Hence the lemma ■ 



bj ~ a- Q,- ~ a" _ , , ( V ~ b, )(a f - a,') ; , v/ _ , k 

b i ~ a t b' - ci i (b t - ci i )(b' - cl ) 



From Fig. 1 as well as lemma 1, it is clear that on extended part of membership and non-membership function 
(as represented by dotted line on left side of ai) sum of membership and non-membership values always exceed 
unity. 

And classical definition of membership functions of IF set restricts upper limit at unity. But if constraint 
a < 1 be added to the model, it does not allow objective function Z- ( X ), i = 1 . . . k, to attain values less than a t . 

But since IF goals and/or tolerances are imprecise in nature and objective function z ; (x),i = l...k, is of 
minimization type, DM will be happier to have optimal value of Z- (x), i = 1 ... k, as less than a t as possible. So 
constraint a < 1 may result in less preferable values as optimal solutions to each IF objective function, as shown 
by numerical example 5.1.2. in this paper. Hence it is suggested not to add constraint OC < 1 so that objective 
functions Z- (x), i = 1 . . . k, may attain values as less than a t as possible. 

On the other hand, the constraint (X > (5 suggests that minimal degree of acceptance cannot be less than 
maximum degree of rejection. This constraint further implies that objective function z t (x) may not attain 
values greater than c ( C ; being abscissa of the point of intersection of membership function and non- 
membership function of objective function z, (x), i = 1 . . . k ). But since IF goals and/or tolerances are imprecise 
in nature, in many cases such constraint may make the problem infeasible, as shown by numerical example 
5.2.1. Hence it is suggested to remove the constraint (X > /3 so that optimal solution may exist with objective 
function z,(x),i = l...k, attaining values greater than C ; . 

And constraint OC > 0 does not allow any objective function z,(x),i = l...k, to take value greater than b t . 
But since IF goals and/or tolerances are imprecise in nature, in many cases such constraint may make the 
problem infeasible, as shown by numerical example 5.2.2. in this paper. Hence it is suggested to remove 
constraint OC > 0 so that optimal solution may exist and objective function(s) Z ; (x),i = l...k, may attain 
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value(s) greater thanZz . It is analogous to removal of upper limit of membership function of fuzzy set theory as 

suggested by Wu et al (2015). Therefore classical membership function now transforms into one strictly 
monotonic function, which seems more useful in decision making under IF environment. 

And classical definition of non-membership functions of IF set restricts upper limit at unity. If such constraint 
/?<lbe added to the problem, it may lead to an infeasible solution, as shown by numerical example 5.2.3. 
Hence it is suggested not to add constraint f) < 1 . Similarly removal of constraint p > 0 is suggested. Therefore 
classical non-membership function now transforms into one strictly monotonic function, which seems more 
useful in decision making under IF environment. 

Hence in an MOLPP under IF environment, it is suggested to use strictly monotonic T (+) -characteristic 
functions in place of membership functions and strictly monotonic T^’-characteristic functions in place of non- 
membership functions; and hence model (3) may be rewritten as follows 

max a- p (5) 

subject to 

T™ (zpx)) > a,i = l...k, 

7j. w (z.(v))</?,i= l...k, 
x s X,a, P unrestricted in sign. 

The task is to find an I-Pareto-optimal solution to model (5) and hence Pareto-optimal solution to MOLPP (1) 
under IF environment. To find I-Pareto-optimal solution to problem (5), the following lemmas are useful. 

Lemma 2 If x* G X is unique optimal solution to model (5), then x* is Pareto-optimal solution to the MOLPP 
(1) under IF environment. 

Proof Let x* G X be unique I-Pareto-optimal solution to model (5). If possible let x * be not Pareto-optimal 
solution to model (1) under IF environment. Then from definition of Pareto-optimality, there exist at least one 

ye. X such that Z ; (v)) < Z,(x*)for all i, i f j and z- (y) < z(x*) for at least one j. Since each T <+) - 
characteristic function T. (+) (z ; ) is strictly monotonic decreasing for minimization type of objective function 
7. V/ — \..k and each ^"’-characteristic function?^ 1 1 (z- ) is strictly monotonic increasing function for 
minimizing type of objective function z ; , i = l...k, we have T l+> (zfx*)) <T. +) (z(y)) V i = 1 . . . k, i ^ j and 
^■ <+) (.Zj (•**)) < Tj (+) (Zj 00 ) and 7^”’ (z i (a*)) > 7^”’ (y)) Vi =l...k,i * j and 7} (-) (z y . (.**)) > Tj~\z j(y)). 

It contradicts that x* is unique optimal solution to model (5). Therefore the assumption is wrong. Hence x* is 
Pareto-optimal solution to MOLPP (1) under IF environment* 

Lemma 3 If x* G X is a Pareto-optimal solution of the MOLPP (1) under IF environment, then x* is an I- 
Pareto- optimal solution of the problem (5) for some < 7j’ + ’ ,7j ( * >, i=l...k. 

Proof Let x* G X be Pareto-optimal solution to MOLPP (1) under IF environment. If possible let x* be not an 
I- Pareto optimum solution to model (5) for any pair of T <+ ’- C haracteristic function and T ( ) -characteristic function 

< 7j. l+) ,7j (_l >,i= l...kto objective function z ; , i = l...k. Then for given functions 7j (+ ’,7j (_) , i = l...k, there 

exists xeX such that 7k <+) (z ; - (x)) > T. <+> (Zj(x*)) for some j and 

Tt + \ Zi {x))>T^\ Zi {x*)),i = \...k,i± ;andr. ( - ) (z ; (x))<r. ( - ) (z.(x*))for some j and 
T.~ ] (z ; (x)) < 7j (_ '(z,(x*)), / = 1. . ,k,i^ j. Since 7J <+) is strictly decreasing function and Tf ’is strictly 
increasing function for each i = l...k, we conclude Z (a) < Z (x*) for some j and z (x) < Z-(x*), 

J J l l 

i = 1 ...k,i V j . It contradicts that x* is Pareto-optimal solution to the MOLPP (1) under IF environment. 
Therefore the assumption is wrong. Hence x* is an optimal solution to model (5) for pair of functions 

< 71 (+) ,7’ (_) > , i = l...k* 
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From lemma 2 and lemma 3, it is clear that if the uniqueness of the optimal solution x* to model (5) is not 
guaranteed, it feels necessary to perform Pareto-optimality test for x*. 

And Pareto-optimality test for x* can be performed by solving optimization problem with decision variables 

x = (x l ,x 2 ,...,x n ) T , Q. = (Q i ,Gl 1 ,...,{\) T andZ5 = (Z5 l ,Z5 2 ,...,Z5 k ) T as follows 

Jc k 

( 6 ) 

1=1 1=1 

subject to 

T i m (z i (x))-Q. i >T i m (z i (x*)),i = l...k, 
r\z i {x))+i5 i <r\z i (x*)),i=\...k, 
xe X,Q,Z5 i >0,i = l...k. 

Lemma 4 Let x,£l and 13 are optimal solutions to model (6). Then 

(1) IfQ = 0 = Z5 j ,\/i = \...k, then x* is Pareto-optimal solution of the MOLPP (1) under IF 
environment. 

(2) If at least one > 0 orX3 i > 0, then x* is not Pareto-optimal solution to MOLPP under IF 
environment (1). Instead of x* X is Pareto-optimal solution to the MOLPP under IF environment (1). 

Proof of analogous lemma is recorded in literature. Also other Pareto-optimality tests exist and may be used in 
place of proposed model (6). 

4 General method to solve multi-objective programming problem under IF environment 

The above ideas can be further integrated into a general framework and an algorithm may be developed to 
obtain I-Pareto-optimal solution and hence Pareto-optimal solution to MOLPP (1) under IF environment. The 
steps of the proposed algorithm are synthesized as follows 

Step 1 . Obtain suitable T (+) -characteristic function for each IF objective in such a way that higher value of T <+> - 
characteristic function always gives better result for objective function. And simultaneously obtain suitable T 1 ' 1 - 
characteristic function for each IF objective in such a way that lower value of T^’-characteristic function gives 
better result for objective function. Well-defined T (+) -characteristic functions and T () -characteristic functions are 
always finite. 

Step 2. Construct single objective optimization problem as model (5). 

Step 3. Solve model (5) and let x* be I-Pareto-optimal solution. 

Step 4. To test Pareto-optimality to this I-Pareto-optimal solution, solve model (6). 

Step 5. Let X, Q. and £5 are optimal solutions of model (3) in step 4. Then two cases may arise: 

(1) If^=0 = U ; — then x* is Pareto-optimal solution of the MOLPP (1) under IF 

environment. 

(2) If at least one Q. >0 orZ3 j >0, then x* is not Pareto-optimal solution to MOLPP (1) under IF 

environment. Instead of x*, x is Pareto-optimal solution to the MOLPP (1) under IF environment. 

Step 6. This solution is I-Pareto-optimal as well as Pareto-optimal solution. The algorithm is complete. 

5 Numerical examples 

5.1. Numerical Example 1 



Consider MOLPP under IF environment as follows 
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Table 1 Individual maximum and 


Objective function 


Individual maximum 


Individual minimum 


minimum of objective functions 


Zj (x) = 5-Xj + 5x2 


10 


0 




z 2 (*) - - 8 . 2*2 


3.33 


-14.06 



Table 2 Goals and tolerances 
of objective functions 



Objective 


Goals 


Tolerances for 


functions 


Membership functions non-membership functions 


ZjW 


8 


1.5 2 


z 2 O) 


-2 


2 2.5 



imprecise max Zy (x) = 5xj +5^ (7) 

imprecise max z 2 (x) = 3X| - 8 . 2 x 2 
subject to 

5xj + 7x2 < 12, 9xj + X 2 < 10, — 5X| + 3x2 < 3, Xj , X 2 S 0. 



To construct classical membership and non-membership functions for each IF objective function zAx),i = 1,2, 
individual maximum and minimum of each IF objective function are computed and given in Table 1. 

5.1.1. Better result on removal of constraint a + J3 <\ 

Now suppose DM specifies goals, tolerances of IF objective functions as in Table 2. Based on these assigned 
goals and tolerances, membership and non-membership functions of each IF objective function may be 
formulated as follows 



A, (^C*)) : 



a(^w)= 



1.5 

0 



3x^ -8.2x2 +4 



, if 5.Xj + 5x2 S 8 

if 6.5 < 5xj +5x 
if 5xj + 5x2 < 6.5 



1 



5x.+5x 0 -6.5 / \ 

1 2 , if 6.5<5xj+5x 2 < 8 andv 1 (z 1 (x)J = 



1 , (/'5-Xj +5 a - 2 < 6 

8 - ( 5x, + 5 x 9 ) 

— , if 6 < 5.Vj + 5x2 ^ 8 

0 , if 5X| + 5x2 S 8 



, if 3.Vj — 8 . 2 x 2 S —2 

, if- 4 < 3 xj - 8 . 2 x 2 < -2 and v 2 ( z 2 (x)) = 
, if 3xj -8.2x2 < -4 



1 



, ;/3xj -8.2x2 < -4.5 



-2-^3xj -8.2x 9 j 



2.5 

0 



if -4.5 < 3x^ -8.2 x'2 S -2 
if 3x^ - 8 . 2 .X 2 >-2 



Using these membership and non-membership functions, MOLPP (7) under IF environment may be converted 
into single objective optimization model as follows [see, e.g. Angelov 1997] 

max a - f5 
subject to 

5x, + 5x 9 - 6.5 3x, -8.2x 9 +4 8-(5x, +5x ? ) -2-(3x, -8.2 x 9 ) 

1 2 > a, 1 2 > a , — v 1 —</3, v 1 —</3, 

1.5 2 2 2.5 

5xj + 7x2 5 12, 9xj + X 2 < 10, - 5xj + 3x2 <3, a > [5, a + J3 < 1, /?, Xp x 2 - 0. 



Using Lingo 15.0.32, 1-Pareto optimal solutions are obtained as follows 
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Table 3 Goals and tolerances 
of objective functions 



Objective 


Goals 


Tolerances for 


functions 




Membership functions non-membership functions 


ZiU) 


8 


1.5 2 


z 2 (x) 


-2 


2 2.5 



a = 1, P = 0, X| * = 0.993, x 2 * = 0.607, z^x*) = 8, z 2 (x*) = -2, 

//jUjU*)) = // 2 (z 2 (x*)) = 1, = v 2 (z 2 (x*)) = 0. 

As discussed in Sect. 3 of this paper, if the constraint a + P < 1 be removed, using Lingo 15.0.32, I-Pareto 
optimal solutions are obtained as follows 

a = 1.131, P = 0, x { ** = 1.045, x 2 ** = 0.594, z^x**) = 8.197, z 2 (x**) = -1.738, 

^jfzjlx**)) = 1.131,7/ 2 (z 2 (x**)) = 1 . 13 1. Vj (zj (x * *)) = -0.098, v 2 (z 2 (x * *)) = -0.105, 

To test Pareto optimality of these I-Pareto optimal solutions, the following problem is now solved 

2 2 

max Z h + I (7 ■ 

(=1 (=1 

subject to 

5x, + 5x 9 — 6.5 3x, - 8.2x 9 + 4 

— 1 z ^ > 1.1311476, — — > 1.1311476, 

-2 - (3x, - 8.2x 9 ) 

+ U, < -0.0983608, — — + U 9 < -0.1049180, 

2 1 2.5 z 

5.Vj + 7 x 9 < 12, 9X| + x 9 < 10, - 5x^ + 3x 9 < 3, 

X| , x 2 >0, Vi=l,2. 

Using Lingo 15.0.32, Pareto optimal solutions x , £1 and 13 are obtained as follows 

D, = 0, = 0 = U 1 = U 2 , X[ =1.045, x, =0.594, Zl (x) = 8.197, z 2 (x) = -1.738 
// 1 (z 1 (x))= ji 2 (z 2 (x))= 1.131, ^(z^x)) = -0.098, k 2 ( z 2 ( x ))= -0.105. 

Therefore it is observed that removal of constraint a + P < 1 generates more acceptable I-Pareto optimal 
solutions with both objective functions attaining better values having higher degrees of acceptances as well as 
lower degrees of rejections. Hence these solutions X** , z(x**)i.e. X, z(x)are more acceptable to DM and 
they have the additional property of being Pareto optimal. 

5. 1.2. Better result without adding constraint a < 1 (so constraint a + P < 1 is not present) 

Again suppose that DM specifies goals and tolerances for each IF objective function as given in Table 3. Based 
on these assigned goals and tolerances, membership and non-membership functions of each of the imprecise 
objective functions may be formulated as follows 



1.5 

8 - (5xj + 5x 0 j 






1 , if 5x^+5.^ >8 

5x, +5x, -6.5 

' — - — j-j , if 6.5 < 5xj +5 x 9 < 8 

0 , if 5X| +5x 9 <6.5 
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nU (■*))= 



8-(5.Vj +5xn 



0 



, if5xy + 5 . v 9 <6 



, if 6<5x^ +5 a 9 <1 
if 5aj +5^2 >8 



&(z 2 (x)) = 



1 , if 3xj - 8 . 2 a 2 > -2 

3xj -8.2x2 +4 



2 

0 



if— 4 < 3xj -8.2x-, < -2 
if 3xj -8.2x2 < -4 



V 2 (z 2 (x)) = 



, if 3xy -8.2x2 < -4.5 



- 2 -| 3 .Xj -8.2X2! 



2.5 

0 



if - 4.5 < 3X| - 8.2a'2 S —2 
if 3x^ - 8 . 2 .X 2 >-2 



Using these membership and non-membership functions, MOLPP under IF environment may be converted into 
single objective optimization model as follows [see, e.g. Angelov 1997] 

max a- p 
subject to 

5x, + 5x~ - 6.5 3x, - 8.2x 0 + 4 

— > a, — > a, 

1.5 2 

8-(5x 1+ 5x 2 ) ^ -2-(3 Xi -S.2x 2 ) ^ 

2 ’ 2.5 

5x^ +7x2 <12, 9x^ + X 2 < 10, - 5xj + 3x 2 < 3, 

a > fi, a < 1, Xj , X 2 ^ 0. 

Using Lingo 15.0.32, 1-Pareto optimal solutions are obtained as follows 

a = 1, fi - 0, x\ = 0.993, x* 2 = 0.607, z x (x*) = 8, z 2 (x*) = -2, 

ACziU*)) = M2 (z 2 (x*)) = 1 , VjCzjU*)) = v 2 (z 2 (x*)) = 0. 



As discussed in Sect. 3 of this paper, if constraint OC < 1 be removed, using Lingo 15.0.32, I-Pareto optimal 
solutions are obtained as follows 



a = 1.131, fi = 0, Xj** = 1.045, x 2 ** = 0.594, z x (x**) = 8.197, z 2 (x**) = -1.738, 

// x (z x (x * *)) = 1.131, ju 2 (z 2 (x**)) = 1.131, Vj(zj(x**)) = -0.098, v 2 (z 2 (x**)) = -0.105. 

To test Pareto optimality of these I-Pareto optimal solutions, following problem is solved 
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Table 4 Individual maximum and 
minimum of objective functions 



Objective function 


Individual maximum 


Individual minimum 


zj(x) = 5xj +5 x 2 


10 


0 


z 2 (x) = 5xj +x 2 


6 


0 


zj (x) = 3xj - 8.2x 2 


3.33 


-14.06 



Table 5 Goals and tolerances 
of objective functions 



Objective 

functions 


Goals 




Tolerances for 


Membership functions non-membership functions 


Zi(x) 


7 


1.5 


2 


z 2 (x) 


2 


2 


2.5 


z 3 (x) 


-2 


2 


2.5 



2 2 

max Z ^ + Z O- 
t=l i = 1 

subject to 

5x, + 5x 9 -6.5 3x, — 8.2x 9 +4 

— a >1.1311476, — O, > 1.1311476, 

1.5 1 2 z 

8 - (5x, + 5 Xr, ) -2 - ( 3x, - 8.2x 9 ) 

— — + U, <-0.0983608, — — + 0 9 <-0.1049180, 

2 1 2.5 Z 

5xj + 7 .X 2 ^ 12, 9x^ + x 9 <10, — 5X| + 3 .X 2 ^ 3, Xj , X 2 , Qj , 15 > 0, Vi=l,2. 



Using Lingo 15.0.32, Pareto optimal solutions x , Cl and 15 are obtained as follows 



C\ =£2, = 0 = 0; = U 2 , x, =1.045, x 2 = 0.594, z^ (x) = 8. 197, z 2 (x) = - 1.738 
//j (zi (x)) = ju 2 (z 2 (x)) - 1.131, P 1 (z 1 (x)) = -0.098, k 2 (z 2 (x)) = -0.105. 

Therefore it is observed that removal of constraint OC < 1 generates more acceptable I-Pareto optimal solution to 
MOLPP (7) under IF environment in which both objective functions have better values having higher degrees of 
acceptances as well as lower degrees of rejections. Hence these solutions X** , z(x**) i.e. X, z(x) are more 
acceptable to DM and they have the additional property of being Pareto optimal. 

5.2. Numerical Example 2 

Consider MOLPP under IF environment as follows 



imprecise max 5X| + Sxj, imprecise min 5xj + Xj, imprecise max 3x^-8.2x 9 (8) 

subject to 5xj + 5x^ < 12, 9x^ + X 2 < 10, -5xj + 3x2 < 3, Xj , X 2 > 0. 

To construct classical membership and non-membership functions for each IF objective function 
Zj (x), i = 1, 2, 3 , individual maximum and minimum of each IF objective function are computed and given in 
Table 4. 

5.2.1. Solution found on removal of constraint a > P 

Now suppose that DM specifies goals, tolerances of imprecise objective functions as in Table 5. 

Based on these assigned goals and tolerances, membership and non-membership functions of each of the IF 
objective functions are formulated as follows 
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A 1 (ziW) = i 



5->fj + 5x 0 -5.5 



1.5 

0 



, if 5.Tj + 5x2 ^ 7 

, if 5.5 < 5xj +5 x 2 < 7 and^M) = j 
, if 5xj +5 x 9 <5.5 



1 



, !/ 5 xj + 5 x 2 



<5 



7-(5xj +5x 9 j 



ji 2 {z 2 {xj) = \ 



1 , if 5xj + x 2 < 2 

4 -( 5 xi + x 2 ) 

2 



and v 2 (z 2 (x)) = j 



ff(z 3 (x)) - 



0 
1 

3X| -8.2x' 2 +4 



2 

0 



if 2 < 5xj + x 2 < 4 
, if 5x^ + x 2 > 4 
i / 3 X | - 8 . 2 x 2 > -2 

if- 4 < 3xj - 8.2x 2 < -2 and v 3 (z 3 (x)) = <j 
z/ 3xj -8.2 x 2 < -4 



if 5 < 5x^ + 5 x 2 < 7 
0 , if 5X| + 5x 2 > 7 

1 , ;/5X| + x 2 > 4.5 



(5x 1+ x 2 )-2 



2.5 

0 

1 



-2-|3xj -8.2x 2 j 



2.5 

0 



if 2 < 5X| + x 2 < 4.5 
, if 5xj + a - 2 < 2 
, z/3xj -8.2 x 2 < -4.5 

if -4.5 < 3X| -8.2x 2 < -2 
, if 3xj - 8 . 2 x 2 >-2 



Using these membership and non-membership functions, MOLPP (8) under IF environment may be converted 
into one single objective optimization problem as follows [see, e.g. Angelov 1997] 



max a — ft 
subject to 



5xj +5 x 2 -5.5 



1.5 



> a. 



-(5xi+x 2 ) 



3x, - 8.2x ? + 4 
>a , — >a , 



7-(5x 1+ 5x 2 ]^ (5x, +x 2 )- 



2 -2 - (3x, — 8.2x 0 ) 

— < /], _ v 1 L ’ 



^ A 



2 2.5 2.5 

5xj +7 x 2 < 12, 9xj +x 9 < 10, -5xj +3x 2 <3,a> j3, a+/3 < l,A x p x 2 > 0. 

Using Lingo 15.0.32, it is found that the problem has no feasible solution. As discussed in Sect. 3 of this paper, 
if the constraint that a > fi be removed, using Lingo 15.0.32, I-Pareto optimal solutions may be obtained as 
follows 

a = 0.284, p = 0.572, Xj = 0.561, x 2 = 0.624, z, (x*) = 5.927, z 2 (x*) = 3.431, z 3 (x*) = -3.431, 

M\(Z\ (**)) = // 2 (z = = 0.284, Kjlzjfx*)) = 0.537, v 2 (z 2 (x*)) = v^iz^ix*)) = 0.572. 

To test the Pareto optimality of these I-Pareto optimal solutions, the following problem is solved 



3 3 

max Z + Z U 

z=l z'=l 

subject to 



■ ^ > 0.2844037, 



+ Uj <0.5366972, 



4 ~ (Sxj + x 2 ) 

2 

^5xj + x 2 ) - 2 



■O, >0.2844037, 



5xj + 5 x 2 - 5.5 
1.5 

7 _ (5xj + 5x^ ) 

2 1 2.5 

5xj + 7 x 2 <12, 9x| + x 2 < 10, - 5xj + 3x 2 < 3, Xj , x 2 , , 0 ; - >0, Vi=l,2,3. 



3x^ — 8.2 x 9 + 4 



■ > 0.2844037, 



+ U 0 <0.5724771, 



-2-^3xj -8.2 x 2 ) 



+ U 3 <0.5724771, 



2.5 
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Table 6 Goals and tolerances 


Objective 


Goals 




Tolerances for 


of objective functions 


functions 


Membership functions non-membership functions 




Zi(x) 


8 


1.5 


2 




z 2 (x ) 


1.5 


2 


2.5 




Z 3 (x) 


-2 


2 


2.5 



Using Lingo 15.0.32, Pareto optimal solutions X , £1 and £5 are obtained as follows 



q = Q,=Q = 0 = (”) , =U 2 = U 3 , x, =0.561, x, = 0.624, zSx) = 5.927, ^(x) = 3.431, z 3 (x) = -3.431, 
/^(z 1 (x))= /A(z 2 (x))= /^(z 3 (x)) = 0.284, v 1 (z 1 (x)) = 0.537, v 2 (z 2 (x))= v 3 (z 3 (x)) = 0.572. 

Therefore on removal of constraint a> jB , it is found that MOLPP (8) has Pareto optimal solutions under IF 
environment specified by DM whereas presence of the constraint a> jB made problem infeasible. 

5.2.2. Solution found on removal of constraint a > 0(so constraint a > fB is not present) 

Again suppose that DM specifies goals, tolerances of each IF objective function as in Table 6. Based on these 
assigned goals, and tolerances, membership and non-membership functions of each IF objective function are 
formulated as follows 



A(mW) = 



1 , ;/5X| +5.x'2 > 8 

5;^ + 5x 9 -6.5 



1.5 

0 



M 2 {z 2 (x)) = 



- , if 6.5 < 5a'| + 5x 2 < 8 and v, ( z, ( x )) = 
, if 5xj +5x2 < 6.5 
1 , if5x^ +x 2 < 1.5 



if 5x, + 5x~ < 6 



8-^5xj +5 x 9 ) 



3.5 -|5x^ +X 2 j 



/i,(z 3 (x)) = 



3x^ -8.2x2 +4 



if 1 . 5 < 5 x 1 + x 2 < 3 . 5 andv 2 (^W) = 

, if 5xj + x 2 > 3.5 
, i/3x'| -8.2x2 > -2 

, if- 4 < 3xj -8.2x 2 < -2 and V 3 (z,(x)) = 

, if 3X| -8.2x2 <— 4 



if 6 < 5xj +5 x 9 < i 
if 5xj +5 x 9 > 8 
, if5x 1 + x 2 > 4 



(5x^+x 9 j-1.5 



2.5 

0 



1 



if 1 .5 < 5X| + x 2 <4 
, if 5x^ +X 2 <1.5 
, ifix^ -8.2x2 < -4.5 



-2-(3.Yj-8.2x-,) 



2.5 

0 



if -4.5 < 3X| -8.2x2 <-2 
if 3.Xj -8.2x2 >—2 



Using those membership and non-membership functions, MOLPP (8) under IF environment may be converted 
into single objective optimization problem as follows [see, e.g. Angelov 1997] 



max 

subject to 

5xj +5 x'2 
1.5 

8 ~(5xj + 
2 

5xj + 7 x'2 



a- fj 



6.5 3.5-(5xj+x 2 ) 3xj-8.2x 2 +4 



> a , 



> a. 



> a , 



' >x 2 ) „ „ ( 5x l + * 2 ) " L5 „ 0 ~ 2 - ( 3x l " 8 - 2 x 2 ) 



<p. 



*P, 



-ZP, 



2.5 2.5 

< 12, 9xj + Xj < 10, — 5xj + 3.x-, < 3, cc > p,a+p < 1, P, Xj, x '2 S 0. 
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Table 7 Goals and tolerances 
of objective functions 


Objective 

functions 


Goals 


Tolerances for 

Membership functions non-membership functions 




Zi(x) 


9.5 


1.5 


2 




z 2 (x) 


1.5 


2 


2.5 




z 3 (x) 


-2 


2 


2.5 



Using Lingo 15.0.32, the problem has no feasible solution. If constraint OC > Obe removed, problem becomes 

max a- p 
subject to 

5x, + 5x 0 - 6.5 3.5 - ( 5 a, + x~ ) 3x, - 8.2x 9 + 4 

— 1 > a, — — >a , — 1 >a, 

1.5 2 2 

8-(5x,+5x 0 ) (5x, + x~)-1.5 —2 — (Sx, - 8.2x 0 ) 

- I ±L < p, - i ZZ < ^ - i ±L < ft' 

2 2.5 2.5 

5xj +7x2 <12, 9xj + x 9 <10, - 5xj + 3x2 < 3, a+P < 1, P, Xj, X2 > 0, a unrestricted in sign. 

Using Lingo 15.0.32, 1-Pareto optimal solutions are obtained as follows 

a = -0.046, P = 0.837, x,* = 0.576, x 2 * = 0.710, z^x*) = 6.43, z 2 (x*) = 3.59, z 3 (x*) = -4.094, 
/tjlZjlx*)) = /i 2 (z 2 (x*)) = // 3 (z 3 (x*)) = -0.046, VjfZjCx*)) = 0.785, v 2 (z 2 (x*)) = ^3(23 (•**)) = °- 837 - 

To test the Pareto optimality of these I-Pareto optimal solutions, following problem is solved 
3 3 

max E !!■+ I ?5 • 
i=l ' i=l 1 

subject to 

5x, + 5x 0 - 6.5 8-1 5x, + 5x 9 ) 

— i Z -ft. >-0.0458716, — — + U , < 0.7844037, 

1.5 1 2 1 

3.5 - ( 5x. + x 9 ) I5X|+x 9 1 — 1.5 

— — - ft, > -0.0458716, - — — + U, < 0.8366972, 

2 2.5 

3x, — 8.2x, +4 — 2 — (3x. -8.2x^) 

— + -ft, > -0.0458716, — — + G, < 0.8366972, 

2 i 2.5 J 

5xj+7x 2 <12, 9x|+x 2 <10, — 5x^ + 3 x 2 < 3, x^ , x 2 , fty , I5 ; - > 0 Vi=l,2, 3. 

Using Lingo 15.0.32, Pareto optimal solutions X , Cl and 15 are obtained as follows 

^ =ft, = 0 = U 1 =U 2 =U 3 , xj =0.576, T 2 = 0.710, ZjCf) = 6.43, z 2 (x) = 3.59, z 3 (x) = -4.094, 
A 1 (7 1 (^))= M 1 (z 2 (xj)= /t 3 (z 3 (x)) = -0.046, Vjlzjfx)) = 0.785, v 7 (z 2 (x))= v 3 (z 3 (x)) = 0.837. 

Therefore after removal of constraint (X > 0 , it is found that MOLPP (8) has Pareto optimal solutions under IF 
environment specified by DM whereas presence of the constraint Ot > Omade problem infeasible. 

5.2.3. Solution found without adding constraint p < 1 (so constraints a> P and a > 0 are not present) 

Again suppose that DM specifies goals, tolerances of IF objective functions as in Table 7. 
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Based on these assigned goals and tolerances, membership and non-membership functions of each of the IF 
objectives are formulated as follows 



A(ziW) = j 



5a'j +5x-, -8 



M z 2 (■*)) = 



, if5xy + 5x 9 > 9.5 

, if 8 < 5X| + 5x 2 < 9.5 - U ( m W) = 

, if 5.Vj + 5^2 < 8 
, if 5xj + x 9 <1.5 

(5.^ +x 9 ) ^ ^ + ^ ^ ; v 2 (z 2 (x)) = { 



1.5 

0 

1 



1 ,if5x-y+ 5x 2 <7.5 
9.5 — (5xfj + 5x 9 ) 



2 

0 



if 7.5 < 5X| + 5x 9 <9.5 
, if 5xj + 5x 9 > 9.5 
if5x^ + x 9 > 4 



( 5 xi +^ 2 )- 1 - 5 



/^(z 3 W) = 



if 1.5 < 5xj +x 2 < 3.5 ’ U 2 W] - y 1 ^ , if 1.5 < 5xj +x 2 < 4 

0 , if 5X| +x 2 >3.5 0 , if 5xj+.v 9 <1.5 

1 ,if3x l - 8.2x 2 >-2 [ 1 ,i/'3x | -8.2x 2 <^1.5 

3x, -8.2x 9 +4 , -2-hx.-S.2xr,) 

— 1 = — — , if- 4<3aj-8.2x 2 <-2. v 3 (z 3 (x))=j LJ 2^ ; j _4 



2 

0 



if 3xj -8.2 x 2 < —4 



2.5 

0 



5<3x^-8.2x 9 <-2 
if 3x^ -8.2x 2 >-2 



Using those membership and non-membership functions, MOLPP (8) under IF environment is converted into 
single objective optimization model as follows [see, e.g. Angelov 1997] 



max a- ft 
subject to 

5x, + 5 x 9 - 8 3.5 — 1 5x, + x 9 ) 3x, - 8.2x 9 + 4 

— >a, — >a, — 1 > a, 

1.5 2 2 



9.5-(5x 1+ 5j 2) (Sxj +x 2 )- 1 -5 ^ -2-(3x 1 -8.2x 2 ) 

2 ~ 2.5 ~~ 2.5 

5xj + 7 x 2 < 12, 9xj + x" 2 < 10, - 5xj + 3x 2 < 3, 

P < 1, a+P <1, p, Xj , x 2 > 0, a unrestricted in sign. 






Using Lingo 15.0.32, it is found that the problem has no feasible solution. 

As discussed in Sect. 3 of this paper, if the constraint p < 1 be not added, the problem becomes 



max a- p 
subject to 

5x, + 5 x 9 - 8 3.5 — 1 5x, + x 9 ) 3x, - 8.2x 9 + 4 

— 1 >a, — —>a, — 1 > a, 

1.5 2 2 

9.5 — ( 5x, + 5x 9 ) ( 5x, + x 9 ) - 1 .5 -2 - ( 3x, - 8.2x 9 ) 

v 1 < p, v 1 ±L < p^ v 1 < p 5 

2 2.5 2.5 

5xj + 7 x 2 < 12, 9xj + x' 2 < 10, - 5xj + 3x 2 < 3, 

a+P <1, p, Xj, x 9 > 0, a unrestricted in sign. 

Using Lingo 15.0.32, 1-Pareto optimal solutions are obtained as follows 
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a = -0.349, P = 1.079, xf* = 0.675, x 2 * = 0.820, z { (x*) = 7.477, z 2 (x*) = 4.197, z 3 (x*) = -4.697. 
/e l (z 1 (x*)) = ie 2 (z 2 (x*)) = ^ 3 (z 3 (x*)) = -0.349, v 1 (z 1 (x*)) = 1.011, v 2 (z 2 (x*)) = v 3 (z 3 (x*)) = 1.079. 

To test the Pareto optimality of these I-Pareto optimal solutions, following problem is now solved 
3 3 

max £ Q- + £ U - 
i = 1 ' i = 1 ' 
subject to 

5x. + 5x 0 -8 9.5 - ( 5x, + 5x 0 ) 

— O, >-0.3486239, — — + U, <1.011468, 

1.5 1 2 1 

3.5 - 1 5x, + x 2 ) (5x,+x 9 1 — 1.5 

— — - a, >-0.3486239, - — — + U- < 1.0788999, 

2 ” 2.5 2 

3x, -8.2x 0 +4 -2 - ( 3x, - 8.2x 0 ) 

1 2 O, >-0.3486239, — — + 0-, <1.0788999, 

2 > 2.5 3 

5X| + 7 % 2 < 1 2, 9 xj+x 2^10 , -5xj + 3+2 < 3,Xp x-, ,£1,U ; . > 0, Vi=l,2,3. 

Using Lingo 15.0.32, Pareto optimal solutions x , £2 and U are obtained as follows 

Q, = a = q, = 0 = U, = U, = U 3 , x, = 0.675, x 2 = 0.820, z, (x) = 7.477, z 2 (x) = 4.197, z 3 (x) = -4.697, 
//i(Zi(x))= // 2 (z 2 (x))= J « 3 (z 3 (x)) = -0.349, K 1 (z 1 (x)) = 1.011, V 2 (z 2 (x)) = V 3 (z 3 (x)) = 1.079. 

Therefore if constraint ft < lbe not added, it is found that MOLPP (8) has Pareto optimal solutions under IF 
environment as specified by DM whereas presence of the constraint ft < 1 may make problem infeasible. 

5.3. Numerical application of proposed algorithm in large scale steel plant 

5.3.1. Usage of T(+)-characteristic functions and T(-)-characteristic functions in steel-iron industry problem 

This example is based on paper titled A multi-objective model for purchasing of bulk raw materials of a large- 
scale integrated steel plant by Zhen Gao, Lixin Tang (2003) published in International Journal of Production 
Economics. Instead of assigning crisp goals of the objective functions, IF goals are assigned to objective 
functions. Moreover in place of membership and non-membership functions, T (+ ’-characteristic functions and T ( " 
’-characteristic functions are used in this MOLPP with IF goals. Here the problem is of purchasing raw materials 
of a large scale steel plant. In steel-iron industry, selection of appropriate items is the key to reduce production 
cost. Decision on quantities of raw materials using strong professional knowledge on steel iron metallurgy is 
another key objective. Selection of vendors to keep stability and quality of supply of raw materials is also key 
objective in steel-iron industry. Using two dimensional vectors x- denoting order quantity of j"‘ item of raw 

materials from i' h vendor, and assuming the model as single time phase model, the constraints under 
consideration may be taken as follows: purchasing budget constraint, production demand constraint, inventory 
capacity constraint, technological requirement constraint, vendor resource constraint. Suppose the problem is to 
have seven vendors and thirteen items that belong to four large kinds of bulk raw materials- ore, lump ore, 
pellet, coal [see Gao 2003], And analogous to Gao (2003), suppose the MOLPP with IF goals is 

Imprecise min z^(x) = 0.1 1 2Xj ^ + 0.1 27+2 ^ + 0.122x^j + 0. 1 15+^ ^ + 0.1 19Xy ^ + 0.0654x^2 + 0.062 lx,-, 2 (9) 

+ 0.0586xZp +0.0602x'g 9 +0.195x22 + 0.185x'22 +0.09521xj^ +0.0975x24 
Imprecise min z 2 (x) = O.lx^ ^ + 0.155x‘2 j + 0.1 7X4 j + 0.12xj ^ + 0.2 x^ + O.lx^ + 0.25X-,-, +0.15x2-, + 0.3x^ 2 

+ 0.15x22 +0. 12x^2 + 0. IX| 4 +0.15x24 
Imprecise min z^(x) = 0.2Xj ^ + 0. \x^ | + O.I5.X4J + O.nx^ | + 0. 1 3Xj | + 0.2Xj 2 + O.lx-,-, + 0.1 5x^ + 0.22x^ 2 

+ 0. 15.V22 + 0. 1 7x^2 + 0.2x^ + 0.15+24 
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Table 8 Purchasing decision of bulk raw materials of a large-scale integrated steel plant with IF goals 



Cases 


I 






II 






III 






Objective functions 


ZiU) 


z 2 (x) 


z 3 (x) 


ZiU) 


Zo(x) 


z 3 (x) 


ZiO) 


z 2 (x) 


z 3 (x) 


Goals 


16 


18.5 


27.5 


15.5 


22.5 


23.8 


15.1 


20.4 


24.4 


j/ + W)) 


0.3 


0.3 


0.3 


0.5 


0.5 


0.5 


0.4 


0.4 


0.4 


Tolerances 




















of 3 


0.5 


0.5 


0.5 


0.7 


0.7 


0.7 


0.6 


0.6 


0.6 


Solutions by 


No feasible solution found 


No feasible solution found 


No feasible solution found 


Angelov’s technique 




















I-Pareto optimal 


15.833 


18.333 


27.333 


15.829 


22.829 


24.129 


15.734 


21.034 
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subject to 
Zj(je) < 16.373, 

1 . 2 xj i + 0 . 9 X 31 +x 4i + 1-1 x^ | +0.95x^j > 60, 

1. 25x^2 +0.95 x2 9 +I.I 5 .X 32 + 1.05xg2 -30. 

I. 3 .T 33 +I.I.X 33 > 10, 1.12xj^ + 1 . 24 X 3 ^ > 70, 

2x^2 + 2^29 + 2X32 + 2x^2 + 3 xgj + 3.X33 = x^ j + X3 ^ 

x.. > 0,Vi = 1...7, / = 1...4. 

y 



+ X41 +x 5 f +X71 +Xf 4 +X34, 



First individual maximum and minimum of each objective function zfx),i = 1,2,3 are computed. Based on 
those values, DM suggests goals and tolerances of each objective function z, i (x),i = 1,2,3 as given in Table 8 . 
Again using these goals, tolerances; strictly monotonic decreasing T (+ ’-characteristic function and strictly 
monotonic increasing T ( ) -characteristic function of each objective function z ; (x),i = 1,2,3 are constructed. 
Different cases may be considered based on different goals and tolerances (assuming that these are supplied by 
DM or presumed values) and the results are tabulated in table 8 . 

On Table 8 , in case I, case II as well as case III, it is shown that there exists no optimal solution if usual 
membership and non-membership functions are used and if existing techniques are followed in optimization 
technique under IF environment. But if strictly monotonic T (+) -characteristic functions and ^‘’-characteristic 
functions are used, proposed algorithm generates I-Pareto optimal solutions and hence Pareto optimal solutions. 
Moreover the Pareto optimal solutions in case I, case II and case III in Table 8 are more acceptable to the DM 
than crisp solution as given by Gao (2003). In proposed algorithm, the optimal values of Zj(x),z 2 (x),z 3 (x) in 
case I, case II and case III are Pareto optimal whereas the optimal values by Gao (2003) are not. And result in 
case I is better than solution number 2, case II is better than solution number 3and solution number 4, case III is 
better than solution number 1 in crisp optimization by Gao (2003). 

5.3.2. Another set o/T <+ ’-characteristic functions and T^-characteristic functions 
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It is discussed in this paper that traditional membership function is not the only function whose higher value 
gives better value of objective function. Similar statement holds for non-membership function as well. Now T <+> - 
characteristic functions and T^-characteristic functions in example 5.3 are constructed analogous to traditional 
membership and non-membership functions. But for minimizing type of objective functions z ; (x) , one example 
of T (+) -characteristic function may be -z ; (x)and of ^'’-characteristic function may be z,-(x) itself. These are 
useful specially when goals and/or tolerances cannot be defined easily by DM or even a DM may not be present 
at all; so that traditional membership and/or non-membership functions cannot be defined for IF objectives. 
Using -z,-(x)as T <+ ’-characteristic function and z,-(x)as ^‘’-characteristic function for each IF objective 
z,- (x) Vi = 1,2,3 and applying proposed algorithm on MOLPP (9) under IF environment, single objective 
optimization problem is obtained as 

max a - ft 
subject to 

—Zj(x) > or, Vi = 1, 2,3, 

Zj(x)< j3,\/i = 1,2,3, 

Zj(x)< 16.373, 

1 . 2 xj i + 0.9x^^ + x^i + l-lx^j + 0.95x.yi > 60, 

I ,25x|2 +0.95x 22 + 1-15x22 + 1 -05x^9 >30, 

1 .3x^2 + I ■ I +53 -10, 1.12xj^+1.24x24 >70, 

2 xp + 2 x 29 + 2 x 22 + 2 xg 9 + 3x 33 + 3x ^ 3 = X| | + x 3 ^ + X 4 1 +x^j +Xr ]\ + x i4 +x 34’ 
x-j > 0, Vi = 1 ... 7, j = 1 ... 4, a, P unrestricted in sign. 

Using Lingo 15.0.32, 1-Pareto optimal solutions are obtained as follows 



a = -23.704. p = 23.704, z,(x*) = 16.373, z 2 (x*) = 23.704, z 3 (x*) = 23.704. 
To test the Pareto optimality of these I-Pareto optimal solutions, following problem is now solved 
3 3 

max X £> + X 75 • 
i=l ' i=l ' 

subject to 

-Zj(x)-D, >-16.373, z 1 (x) + 75 1 < 16.373, 

-z 2 (x)-£2, >-23.70491, z,(x) + 75 2 < 23.70491, 

-z 3 (x)- £2, >-23.70491, z,(x) + 75 3 <23.70491, 

Zj(x) < 16.373, 1.2xj y +0.9x 3 j +x^^ + l.lx 3 i + 0 . 95 x 7 j ^ 60, 

1. 25x^2 +0.95x22 +1-15x29 + 1 - 30, 

I. 3 .V 22 + l-l x 53 — 10’ 1-12 x j4 +1.24x24 — ^0, 

2x^2 + 2 x 29 + ^ x 32 + ^ x 62 + ^ x 33 + ^ x 53 = X 11 +x 31 +x 41 +x 51 +x 71 +x 14 +x 34’ 
x.j > 0 , Vi = 1 ... 7 , 7 = 1 ... 4, £2^,75. > 0 , Vi=l,2,3. 



Using Lingo 15.0.32, Pareto optimal solutions x , £2 and 75 are obtained as follows 



£2, =£2, =£2, =0 = 75j =75, =75 3 ,z 1 (x) = 16.373, z 2 (x) = 23.704, z 3 (x) = 23.704. 

This is another set of Pareto optimal solutions. Moreover instead of these T <+) -characteristic functions and T ( ‘’- 
characteristic functions, different type of functions may be constructed as necessary. 
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6 Conclusion 

Technique to find Pareto optimal solutions to multiple objective linear programming problems under IF 
environment is introduced in this paper. IF technique is one of the richest apparatus for formulation of 
optimization problems under imprecise environment, thereby generating more satisficing result than crisp 
optimization technique. It is shown in this paper that all of the constraints in existing techniques are not 
necessary. Presence of those constraints may make a problem infeasible. These constraints were necessary due 
to definition of membership and non-membership functions of IF goals. Consequently membership and non- 
membership functions of IF goals cannot be used in optimization technique under IF environment. Strictly 
monotonic T (+) -characteristic functions and T ( ) -characteristic functions are proposed to be used in place of 
existing membership and non-membership functions of IF set. 

Moreover, in multiple objective optimizations, Pareto-optimality of solutions guarantees rationality of the 
decision. One new technique for finding Pareto optimal solutions to MOLPP under IF environment is introduced 
and discussed in this paper. Hence the proposed algorithm generates I-Pareto optimal solutions, followed by 
Pareto optimal solutions to MOLPP with IF goals. 

As shown in the application, proposed algorithm generates more acceptable solutions than crisp optimization 
technique by Gao (2003). And this is a general procedure to obtain an I-Pareto-optimal solution that has the 
additional property of being Pareto-optimal. Proposed algorithm can find Pareto optimal solutions to MOLPP 
with IF goals even when goals and/or tolerances cannot be set by using T <+ -characteristic function and T 1 ' 1 - 
characteristic function, as shown in example 5.3.2. 

From this it may further be concluded that, under these circumstances, the issue of getting no solution even by 
using IF set theory yielded from several published methods seems worthwhile or even necessary to reconsider. 
The same is also true for nonlinear problems of this type. 
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